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Creating maximally entangled atomic states in a Bose-Einstein condensate 
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We propose a protocol to create maximally entangled pairs, 
triplets, quartiles, and other clusters of Bose condensed atoms 
starting from a condensate in the Mott insulator state. The 
essential element is to drive single atom Raman transitions 
using laser pulses. Our scheme is simple, efficient, and can be 
readily applied to the recent experimental system as reported 
by Greiner et al. [ Nature 413, 44 (2002)]. 
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The physics of quantum degenerate atomic gases con- 
tinues in its rapid pace of development, and remains 
one of the most active research areas in recent years 
[1]. Increasingly, theoretical and experimental attentions 
are directed towards the underlying quantum correlation 
properties of the condensed atoms. It seems likely that 
such quantum states of matter might prove to be a fertile 
ground for exploring quantum information science appli- 
cations. 

Recently, a quantum phase transition was observed in a 
system of Bose condensed atoms immersed in a periodic 
array of optical potentials [2]. As expected, when ex- 
pressed in the simple Bose- Hubbard form [3] , the ground 
state of such a system is controlled by essentially two 
parameters: 1) the on-site atom-atom interaction u for 
atoms in the same spatial mode of each individual optical 
well; and 2) the nearest neighboring well (single) atom 
tunnelling rate J (taken as positive). When J ^ \u\, the 
condensate ground state is in the usual superfluid (delo- 
calized single atom) state. On the other hand, a Mott 
insulator state arrives in the opposite limit |m| ^ J. In 
a Mott state, atoms are localized inside individual wells. 
The condensate ground state takes the form of a product 
of Fock states with an integer number of atoms on each 
site. The transition from superfluid to Mott insulator is 
predicted to occur at \u\/ J > z x 2.6 with z the number 
of nearest neighbors in the periodic well lattice [3,4]. 

The experimental system that yielded the first clear 
demonstration of the superfluid/Mott-insulator transi- 
tion enables individual tuning of the values for both J 
and u [2]. In the experiment, the average occupations 
per well was around 1-3 atoms, which could potentially 
form elementary building blocks for atomic qubit based 
quantum computing designs [3]. 

In this paper, we propose to create massive maxi- 
mum entangled pairs, triplets, quartiles, and other clus- 
ters of Bose condensed atoms in a Mott insulator state. 
The resulting entanglement, with respect to electronic 



internal states, is stable and long lived. In the exper- 
iment [2] ^^Rb atoms in the magnetic trapping state 
\a) = \F = 2,Mf = —2) were used. Other internal 
states can be trapped in the same experimental setup. In 
the simple model to be presented below, a second inter- 
nal state 1 6) that can be coupled to \a) through atomic 
Raman transitions is assumed [2] (as see earlier JILA 
experiments with ^^Rb states \F — 2,Mp = — 1) and 
\F = 1,Mf = 1) [5]). 

In a Mott state, the system dynamics is rather simple 
as there exists a fixed (small) number of atoms within 
each well. If we use the second quantized operators a{a^) 
and h{W) for atoms in the two internal states, the effective 
Hamiltonian for each well can be expressed as [6] 



H = uJi 



(1) 



The second term denotes the single atom Raman coupling 
due to external laser fields with a (real) effective Rabi 
frequency ri(i) [7]. The angular momentum operators are 
the Schwinger representation in terms of the two boson 
modes 



-{b^'a- a^'b), 



-{b^b-a^a). 



(2) 



In the context of SU(2) coherent states of an atomic en- 
semble, these operators have been used extensively for 
discussing spin squeezing and other properties of multi- 
atom nonclassical states [8-11]. In particular, as was 
studied by Molmer and Sorensen [12], an interaction of 
the type uJj generates a maximum entangled N-GHZ 
state [13] starting from all atoms in state \a) or \b). This 
has led to the recent creation of a 4-ion maximum entan- 
gled state [14]. 

Before we discuss our proposal, we summarize the dy- 
namic generation of a maximum entangled state from the 
uJ^ interaction. For simplicity, we assume N is even. A 
maximum entangled N-GHZ state can be written as [12] 



1 ^ A b^^ 

|GHZ)jv = [e'^'^ + e 



m— — 4f 

[e^0i. +e"^a(_l)f-™]|0), (3) 



where new bosonic operators d/c — (& ± a)/\/2 were in- 
troduced along with its inverse b/a= ((i±c)/\/2- is 
the binomial coefficient. Starting from all atoms in state 
|a), i.e. with = at^|0)/\/1v!. The state at time t 

due to a u interaction alone is 



m)) = 



1 



m— - 

e-'"*"'(-l) 



E 



'|0), 



(4) 



where use has been made of Jx = {d)d — c^c)/2. To 
within an overall phase factor \'iP{t)) = |GHZ)jv at ut = 
{2k + l)7r/2 with the shortest time being r — 7r/(2|u|). 
Similarly, starting from state fe^^jO) will also arrive at a 
N-GHZ when ut = {2k + 1)tt/2 [15]. 

How could interaction (1) be turned into the required 
form? Our key observation is that the single atom 
Raman coupling QJy generates nothing but a rotation 
along the y-axis. Therefore, we can effectively rotate 
the J| term into a term. A similar suggestion was 
made recently by Jaksch et. al. [16] in order to tune the 
overall condensate interaction strength to zero (or SU(2) 
symmetric) . 

We therefore suggest operating in a three step protocol 
in the limit when ^ N\u\: 

1) Apply a n/2 pulse 6'(t') = /J" n{t)dt = 7r/2 (of spin 
1/2). During this stage the nonlinear interaction can be 
neglected (because \rt\ ^ ^|^*|)- 

2) Wait for a time \u\t — Tr/2. 

3) Complete the process by applying a — 7r/2 pulse with 
9{t') = —tt/2 [e.g. by arranging for —il or by 
waiting for a 37r/2 pulse as in 1) ]. 

These three steps generate the following effective evo- 
lution, 



[/(2r' -}-r)? 



(5) 



i.e. is rotated by 7r/2 into J^. From a wide range of 
numerical simulations, we find that N-GHZ states with 
extremely high fidelities are realized when > 50 

for (up to 4 atoms). 

While the above scheme works well, it is inherently 
rather slow. In a two component condensate as assumed, 
we denote the 3 relevant scattering lengths as Uaa, a-ab, 
and af,b, and assume that motional ground state to be 
ipooo{r) = exp[-r2/(4a^)]/(V27ra/i)^/^ of a spherically 
symmetric harmonic trap V{r) — Moj'^r'^ /2, we find 



{a-aa + flbb ^ Sflaf,) 



2Trh^ 



1 



M (20Fa,,)3' 



(6) 



with ah = \fhl2MuH the ground state size. For ^''Rb, 
u is very small as Oaa ~ o,ah ~ o,hh- When ujt ^ 
(27r)30 (kHz) as realized in [2], |m| - (27r)20 (Hz) if 
{o-aa + i&f) — 2aah) is of the order of one Angstrom (A). 
It takes approximately 10 (ms) to realize a GHZ state, 



i.e. in a time significantly shorter than the lifetimes from 
both the two-body dipolar [> 6(s)] and the three-body 
inelastic collision [> 200(ms)] losses with less than five 
atoms in each well [17]. 

Another serious experimental concern is that collisions 
can populate Zeemen states other than |a) or For 
most systems, this depopulation also occurs on the time 
scale of ~ It is therefore important to include 

the full manifold of atomic internal states. To this end, 
we consider a spinor-1 condensate of *^Rb atoms in its 
ground state F — \ manifold as realized in the first all 
optical condensate [18]. If aup denotes the bosonic an- 
nihilation operator of state \F — l,Mi? = +,0,—), the 
ground state Hamiltonian within each well becomes 

H'= u{L^ - 2N) 

~ M(a^a^a+a+ + aLaLa-O- 
+2a\_aQa+ao + 2aLala-ao — 2a\_aLa+a- 
-t-2ajaja+a_ + 2a^alaoao), 
with angular momentum type operators [19-21] 
L+= \/2(a^ao + aja_), L_ 



(7) 



(8) 



ajao 



and number of atoms in the well N = a'_^_ 
oLa-. Although seems SU(2) symmetric, it is not be- 
cause La;, Ly, and Lz are not genuine angular moment op- 
erators (for spin-1 atoms); they do not satisfy the Casimir 
relation ^ N{N+1) [22]. As was shown before [19,22] 
multi-atom internal state correlations continue to arise 
dynamically with H' and the addition of single atom Ra- 
man couplings of the type iflfj,,y{ajj^ai, — ala^)/2. Unfor- 
tunately, we have not been able to solve for the combined 
dynamics analytically even for a small number of atoms. 
It is also not apparent how to numerically investigate 
strategies for creating a N-GHZ state in this case. 




4 6 
t (in unit of 7i/u) 

FIG. 1. Solutions of tm and Q. as given by the cross 
points of the two families of curves {utm) = (2A: + 1)-k 
(for k = 0,1,2,3,4) and ^Jl + ■i{Q./uY{ut^) = 2m7r (for 
m = 1,10,30,50, 100). 

Looking back on the two mode model (1) discussed 
earlier, we realize that, with a constant fi, a state with 



two atoms initially in \b) develops into a 2-GHZ state 
within a time of sa 7r/|w|. Specifically, we find 

Cii{t)= -J^e*t*sinm, 

1 M u ~ 1 u 

^ ^ 2 40 2 

1 7; „ - 1 „ 

Co2W= - +i-^e'2*sinm- -e'^tcosOi, (9) 

2 40 2 



with O — Vw^ + 40^/2 for the coefficients of state vector 
expansion 

\m)=C2o{t)^b^^\Q,0)+Co2{t)^a^^\Q,Q) 

+Cii(i)6^a^|0,0). (10) 

In the above Eq. (9), we have omitted a common phase 
factor e~™*. Clearly, Cxi{t) — occurs at 



20t,„ = ^Jl + A{n/uY{utm) = 2m7r. 



(11) 



IViW) becomes a 2-GHZ state (3) when IC20I = IC02I = 
l/\/2. This occurs at = (2fc-|-l)7r since C2o/o2(im) = 
^1 ± g»ut„/2(_;^)m]/2. When |0/u| > 1 both conditions 
can be satisfied at different values of tm and O as shown 
in Fig. 1. The shortest time for a 2-GHZ is then ~ 7r/|u|. 

Based on this observation, we explored numerically 
the dynamics of the Hamiltonian H = u{L'^ — 2N) + 
«0^y(a]^ay — ala^)/2 assuming a constant 0^,^ and all 
atoms initially in the |-|-) state. As expected, we dis- 
covered that maximally entangled states continue to be 
generated at times ~ '^/\u\ for = 2, 3, 4. 

For = 2, we find that we get a 2-GHZ state 
[a^^ + e'-^aj^lo -)|0>0'0>/2 with either a Raman drive 

0+0 or 0_| a controllable phase shift. The 2-GHZ 

state occurs at times of « {2k + l)7r/|M| (/i — 0) or 
{2k + 1)^1 A\u\ (/i = — ) and also times shifted by a small 
multiples of 7r/|0+^| (when |0+^| ^ in their imme- 
diate neighborhoods. The state fidelities are always very 
high as long as k is not too large. 

For Af = 3, only the O^ drive seems to create a 3- 



GHZ state cx (a^^^-|-e'''aL'')|0,0,0) at times differing from 
« {2k -t- l)7r/4|u| by smah multiples of 7r/|0+_|. Maxi- 
mum correlated atomic ensembles in states |-|-) and |— ) 
were previously predicted to occur due to elastic colli- 
sions for a initial condensate in state |0) [25]. 

For A^ = 4, we find that again only the O^ drive al- 
lows for a simple identification of a 4-GHZ state oc (a^ + 

e'V'^)|0,0,0), which also occurs at w {2k + l)7r/4|u| 

and values shifted by a small multiples of 7r/|0^ | in 

its neighborhood. Thus at t sa 7r/4|ii| atoms in wells 
with N = 2 and 4 are both maximum entangled as il- 
lustrated in Fig. 2. In this simulation, we have used 
0+_ = (27r)30 (kHz) and u = (27r)0.25 (kHz). We note 
that their respective values are not important except that 
they scale inversely with the required time. What seems 



to be important is to assure that |0_| /u\ > 100 for up 

to 4-atoms to achieve a high fidelity maximum entangled 
state. 
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FIG. 2. The two oscillating lines are respectively the prob- 
abilities for all atoms in state |-|-) or | — ). Top panel is for 
A'' = 4, while the bottom one is for N — 2 atoms. The verti- 
cal dot-dashed line is at t — 0.5 (ms). 

In conclusion, we have presented a simple and efficient 
protocol for turning a Mott insulator condensate of *^Rb 
atoms in the ground state F — 1 manifold into a source 
for maximally entangled atomic clusters. Our protocol 
is reliable and accessible with current technologies [2] . It 
produces maximum entangled quantum states of Bose- 
condensed atoms with high fidelity. The only noticeable 
drawback seems to be due to the fact that for ®^Rb atoms, 
u oc (02 — ao), i.e. the difference of scattering lengths 
for the two symmetric channels with total spin and 
2. Nevertheless, inelastic decay processes are essentially 
negligible because all spin states of the atomic ground 
state manifold are included. Furthermore, the N-GHZ 
state (|-f)<3^ -t- e"^\-)'^^)/V2 is stable against elastic 
collisions, which are required to conserve the total Mp, 
i.e. atoms in |-|-) (or |— )) state remain in the same state 
after collisions. Thus the slow dynamics is perhaps not a 
major course of concern. Other atomic species (e.g. F = 
1 manifold of ^"^Na [23]) may provider large values of u. In 
Ref. [24], a quantum logic operation between two atoms 
(one each in two neighboring wells) was proposed that 
uses the much stronger (by two orders of magnitudes) 
interaction cx Qab- Application of this in a Mott state 
(with one atom per well) produces GHZ states on a faster 
time scale, although it requires more complicated internal 
state dependent optical trapping. 

Finally, a condensate in a Mott state contains many 
individual wells with identical number of atoms [2] . This 
makes the experimental detection of the entanglement 
(for atoms within each well) relatively easy. One can 
perform the usual parity-type measurement with Ram- 
sey's oscillatory fields technique [14] (again) by driving 
the single atom Raman transition so quickly that collision 



effects are negligible. All wells with the same number of 
atoms thus contribute to the detected signal. General- 
izations of our protocol to more than 4-atoms and other 
related results will be published elsewhere. 

This work is supported by a grant from NSA, ARDA, 
and DARPA under ARO Contract No. DAAD19-01-1- 
0667, and by a grant from the NSF PHY-01 13831. 



[21] T.-L. Ho and S. K. Yip, Phys. Rev. Lett. 84, 4031 (2000). 
[22] O. E. Miistecaplioglu, M. Zhang, and L. You, Phys. Rev. 

A 66, 033611 (2002). 
[23] D. M. Stemper-Kurn et al., Phys. Rev. Lett. 80, 2027 

(1998). 

[24] D. Jaksch, H.-,J. Bricgcl, ,J. I. Cirac, C. W. Gardiner, and 
P. Zoller, Phys. Rev. Lett. 82, 1975 (1999). 

[25] H. Pu and P. Meystre, Phys. Rev. Lett. 85, 3987 (2000); 
L.-M. Duan et al., ibid. 3991 (2000). 



[1] See the extensive Ust of references at 

http: / / amo.phy.gasou.edu/bec.html/bibliography.html. 
[2] M. Greiner, O. Mandel, T. EssHnger, T. W. Hansch, and 

I. Bloch, Nature 413, 44 (2002). 
[3] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and 

P. Zoller, Phys. Rev. Lett. 81, 3108 (1998). 
[4] S. Sachdcv, Quantum Phase Transitions, (Cambridge 

Univ. Press, Cambridge, 2001). 
[5] M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. 

Hall, M. J. Holland, J. E. Wihiams, C. E. Wieman, and 

E. A. Cornell, Phys. Rev. Lett. 83, 3358 (1999). 
[6] A. S0renscn, L-M. Duan, J. I. Cirac, and P. Zoller, Nature 

409, 63 (2001). 
[7] C. K. Law, H. T. Ng, and P. T. Leung, Phys. Rev. A 63, 

055601 (2001). 

[8] M. Kitagawa and M. Ueda, Phys. Rev. A 47, 5138 (1993); 

Phys. Rev. Lett. 67, 1852 (1991). 
[9] D. J. Wineland, J. J. Bollinger, W. M. Itano, F. L. Moore, 
and D. J. Heinzen, Phys. Rev. A. 46, R6797 (1992). 
[10] C. Orzel, A. K. Tuchman, M. L. Fenselau, M. Yasuda, 

and M. A. Kasevich, Science 291, 2386 (2001). 
[11] Kristian Hclmerson and L. You, Phys. Rev. Lett. 87, 
170402 (2001). 

[12] K. Molmer and A. Sorensen, Phys. Rev. Lett. 82, 1835 

(1999) . 

[13] D. M. Grccnbcrgcr et al., Am. J. Phys. 58, 1131 (1990). 
[14] C.A. Sackett, D. Kielpinski, B.E. King, C. Langer, V. 

Meyer, C.J. Myatt, M. Rowe, Q.A. Turchette, W.M. 

Itano, D.J. Wineland, and C. Monroe, Nature 404, 256 

(2000) . 

[15] For odd values of N, an additional single atom interaction 
of the form QJx (for a duration of tt/2) will be sufficient 
to create a N-GHZ (see Ref. [12]). Such a Raman tt pulse 
can be implemented by taking n — » in in the QJy term 
ofEq. (1). 

[16] D. Jaksch, J. 1. Cirac, and P. Zoller, Phys. Rev. A 65, 
033625 (2002). 

[17] We used the measured two-body and three-body loss 
rates as reported in E. A. Burt et. al, Phys. Rev. Lett. 
79, 337 (1997). Theoretically predicted rates are usually 
smaller. 

[18] M. Barrett, J. Sauor, and M. S. Chapman, Phys. Rev. 

Lett. 87, 010404 (2001). 
[19] C. K. Law, H. Pu, and N. P. Bigclow, Phys. Rev. Lett. 

81, 5257 (1998). 
[20] M. Koashi and M. Ueda, Phys. Rev. Lett. 84, 1066 

(2000). 



